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Abstract 



For the confining phase of SU(2) Yang-Mills thermodynamics we show that the 
asymptotic series representing the pressure is Borel summable for negative (unphys- 
ical) values of a suitably defined coupling constant. The inverse Borel transform is 
meromorphic except for a branch cut along the positive-real axis. The physical pres- 
sure is precisely nil at vanishing and acquires a small imaginary admixture at small 
temperature the latter indicating violations of thermal equilibrium (turbulences). 



1 Introduction and miniature review 



To understand the (thermo-) dynamics of strongly coupled Yang-Mills theories is of 
importance both in view of an unearthing of their interesting mathematical struc- 
tures and an exploitation of these theories in physical applications. 

As discussed in [I][2], SU(2) and SU(3) Yang-Mills theories occur in three phases. 
The deconfming, high-temperature phase, see also [2J H], possesses a nontrivial 
ground state which is composed of interacting (anti)calorons of unit topological 
charge modulus. To analytically access this highly complex dynamics one is lead to 
perform a spatial coarse-graining leading to an effective theory subject to a maximal 
resolution \<p\. Here denotes an emergent, inert, adjoint, spatially homogeneous, 
and BPS saturated scalar field [U El E]- The loop expansion of thermodynamical 
quantities, carried by coarse-grained, topologically trivial modes of potential reso- 
lution smaller than \<p\, converges very rapidly [6]: There is infrared stability due 
to the emergence of mass (adjoint Higgs mechanism), and the action of collective 
quantum fluctuations swiftly decreases with a growing number of loops. 

For a small window of intermediate temperatures a preconfining phase occurs, 
see also [7j. The ground state of this phase is composed of condensed magnetic 
(anti)monopoles. A spatial coarse-graining uniquely generates a scale of maximal 
resolution given by the modulus of a spatially homogeneous, inert, and BPS satu- 
rated complex scalar field (p. The stable excitations in the preconfining phase are 
noninteracting, massive, dual gauge modes, and thus the loop expansion of thermo- 
dynamical quantities is trivial. (It is represented by a ground-state plus a one-loop 
contribution.) 

The objective of the present work is an investigation of the dependence on tem- 
perature T of thermodynamical quantities in the confining phase of SU(2) Yang-Mills 
theory, for detailed discussions see also [HI IS] , where all gauge modes are infinitely 
heavy, the ground state is a condensate of paired, (massless) magnetic center- vortex 
loops, and the propagating excitations are massless or massive spin- 1/2 fermions. 
These solitons are classified according to their topology arising from the number n 
of selfmtersections. The latter give rise to a naked mass nA, A being the Yang- 
Mills scale (or the mass of the stable soliton with n — 1). Naively, that is, when 
neglecting the excitability of internal degrees of freedom within a given soliton and 
when disregarding the (contact) interactions between solitons the total pressure is 
represented in terms of an asymptotic series in powers of a dimensionless coupling 
A = exp(— A/T). Notice that A is strictly smaller than unity for T < A0. 

The question then arises whether, as a matter of principle, sufficient information 
is contained in the asymptotic series to generate the temperature dependence of the 
physical pressure^). The main result of the present work is to answer this question 

1 At T ~ A a Hagedorn transition occurs. The latter goes with a condensation of the selfintersec- 
tion points within densely packed center-vortex loops into a new ground-state: The (anti)monopole 
condensate of the preconfining phase. 

2 By physical we mean that this quantity takes into account the dressing of and the interactions 
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with yes albeit subject to a surprise: Our result predicts the exact vanishing of the 
physical pressure at T = and the perceptible breakdown of thermodynamical equi- 
librium at a sufficiently large temperature. This effect is manifested in terms of an 
imaginary admixture to the real pressure. Once the temperature dependence of this 
real part is known other thermodynamical quantities can be computed by Legendre 
transformations. Our analysis invokes a combination of Borel-transformation and 
analytical-continuation arguments for both A and the Borel parameter. 

The article is organized as follows. In Sec. [5] we discuss and set up an asymp- 
totic series representing the pressure in the confining phase of an SU(2) Yang-Mills 
theory. Modulo an algebraic-in-n uncertainty this series is explicitly known. Sec. [3] 
processes the asymptotic series. Namely, in Sec. 13.11 we perform a Borel transforma- 
tion which recasts the part of the pressure arising from massive excitations into a 
linear combination of polylogarithms. Some of the analyticity structure of the latter 
is discussed subsequently. In Sec. 13.21 we perform the inverse Borel transformation 
numerically for negative arguments of the polylogarithms and observe that the in- 
verse Borel transform is a meromorphic function except for a branch cut along the 
positive-real axis. Motivated numerical fits to the inverse Borel transforms of the 
first few relevant polylogarithms are carried out in Sec. l3.3l for negative values of the 
argument. Subsequently, the fits are continued to positive values of the real part of 
the argument. It is observed that the real parts of the inverse Borel transforms are 
continuous across the cut while the moduli of the imaginary parts are smaller than 
those of the real parts for sufficiently small values of A. An interpretation of this 
result (violation of thermal equilibrium by turbulences) is given in Sec. 13.41 In Sec.H] 
we summarize our results and point towards applications and future research. 

2 Naive series for the pressure 

The excitations in the confining phase of an SU(2) Yang-Mills theory are generated 
by phase jumps and an increase in modulus of a complex scalar field $ which mea- 
sures the expectation of the 't Hooft loop operator [10]. The latter is a dual order 
parameter for confinement. The phase of $ is given by a line integral of the dual 
(abelian) gauge field along an S 1 of infinite radius measuring, by Stokes' theorem, 
the magnetic flux through its minimal surface M.\. The creation of a magnetic 
center-vortex loop is understood as the process of having an infinitely thin flux line 
and its oppositely directed partner (travelling in from infinity) intersect with the S 1 
which leads to their subsequent piercing the surface M.\. This brings into existence 
a propagating, single center-vortex loop. The required energy needed is provided by 
a potential V($) which is initiated by a decaying monopole condensate at the Hage- 
dorn phase boundar^. The potential V(Q) forces the field $ to change its phase 

between naked solitons. 

3 For |$| < A an increase of $'s modulus leads to a decrease of <&'s energy density. At the two 
minima <f> = ± A energy density and pressure vanish. For more information see [IJ [21 [5] . 
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discontinuously in units of n (units of center flux). Collision^ of single center- vortex 
loops lead to twisting and merging thus creating selfmtersections. This process con- 
verts kinetic energy into mass. In reality the very process of thermalization proceeds 
via instable high-mass excitations: Annihilations of oppositely, Z2 charged intersec- 
tion points locally elevate the energy density of the field $. Subsequently, $ ejects 
this energy by phase jumps and an increase of its modulus thus recreating center- 
vortex loops. In the hypothetic case of absolutely stable excitations at arbitrarily 
large selfintersection number n, which once and forever were thermalized to a given 
temperature by a single relaxation process of the field $ to one of its minima, there 
is no ground-state contribution to the pressure. Here we are concerned with this 
hypothetical situation which allows for drastic simplifications compared to the direct 
analytical treatment of the full physical situation and, as we will see below, captures 
sufficient information to recover the latter. 

The total pressure at temperature T is then represented by the following (asymp- 
totic) series (sum over partial pressures of fermionic particles stemming from the 
sectors with n selfmtersections) 

00 00 poo 2 

P = J2 P " = T Y, M " dp log (l + e-^) , (1) 

n=0 ra=0 ^° 71 

where /3 = ^, uj n = \/p 2 + (nA) 2 , A the Yang-Mills scale, and the number M n 
denotes the multiplicity for solitons with n selfmtersections (and naked mass nA): 

M n = 2N n C n . (2) 

Here the factors of 2, of N n , and of C n stand for the spin multiplicity, the number 
of distinct topologies, and the Z 2 charge multiplicity, respectively. One has (two 
possibilities at each intersection point) 

C n = 2 n . (3) 

In Fig.[TJ we list soliton topologies for n < 3. It is obvious that N n represents 
the number of connected vacuum bubble diagrams with n vertices in a A0 4 -theory. 
Modulo algebraic-in-n factors and for large n the number of such diagrams in a 
A0 2-M -theory was found to be 

HM-im^j} Mn (M-ir (4) 

by Bender and Wu [IT] . Thus for Ai = 2 we have 

iV n ~n!!6 n , (5) 



4 Only contact interactions occur due to the complete decoupling of propagating gauge modes 
in the confining phase [TJ [5] . 
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Figure 1: A list of the soliton topologies for n < 3. The number N n is the number 
of connected vacuum bubbles with n vertices in a A0 4 -theory. 

where the ~ sign signals a dependence modulo a factor of the form J2t=i a i n? j 
(L < oo and a/ integer). The total pressure P is then represented by the following 
asymptotic series: 

, 4 / 7^4 



\ n>l • /u 



-/3\/ n' 2 +x 2 



„2 .-.-^v 7 ^" 



n>l 



n>l 



n>l 
L 



* m + ^ S a < B 32A ) n - ! ) • ( 6 ) 

\ Z=0 n>l / 

In Eq. (EI) we have defined /3 = A = e _/3 , and denotes a modified Bessel 
function. From the A dependence of the sought-after physical pressure it should 
later become clear that the very concept of thermalization ceases to be useful for 
T y A or for A f e _1 due to the vicinity to the Hagedorn transition pQ [2J. In 
Eq. ([S]) the first < sign holds strictly for the linear truncation of the expansion of 
the logarithm about unity, and the ~ sign indicates that terms of order (fin)^ 1 
have been neglected in the expansion of the nonexponential factor in the Bessel 
function. The second < sign holds because we have used the large-n expression for 
N n of Eq. ([5]). The coefficients ai, which determine the algebraic factor in N n , are 
unknown at present. 
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Obviously, the sum in the last line of Eq. ([6]) diverges but exhibits the charac- 
teristics of an asymptotic expansion. Notice the formal similarity of the expansion 
in Eq. flB]) with the perturbative loop expansion of the ground-state energy in a 
A0 4 -theory in one dimensional for which Borel summability was proven, see [H] and 
references therein. 



3 Borel transform and physical pressure 

Our strategy to elute the physical pressure from the asymptotic expansion in Eq. ([6]) 
is to perform a Borel transformation of this series. Subsequently, we investigate the 
region of analyticity of the Borel transform. Only for A < is the inverse Borel 
transform real-analytic in A. However, we can analytically continue our results for 
A < to A = Re A ± iO , (Re A > 0) , thus obtaining a prediction of the theory for 
the (perceptible) physical real part of the pressure at physical coupling. 



3.1 Borel transformation 

The Borel transformation removes the factors of n\ in the coefficients of the power 
series in Eq. (jSJ: 

L L n 

P mass (A) = $>£A^!nf+< *3 B Pmass (\) = J2 a iJ2~ Xnnl+l > ( ? ) 

1=0 n >l 1=0 n>l 

where A = 32 A. A sum over n as in Eq. (J7J) defines the polylogarithm Li s (z) for 
complex numbers s and z with \z\ < 1 (here s — — (| +lj and z — A): 



Li -d + oW = E^ l+z - ( 8 ) 



n>l 



By analytical continuation the function Li_/3 + ^(A) is defined for a much larger 

range in A than the definition in Eq. (jBJ) seems to suggest. In any case, there is a 
branch cut for positive-real values of A with A > 1. One has 



.(log(A))-(I + 



ImLi_ r|+0 (A±zQ) = ±tt Y ( _ ( 3 + a } . ( A >!)' ( 9 ) 



> 2 



where T(s) = J °° dt e~* t s_1 the gamma function. Notice that T(— (§ + £)) is finite 
and real due to the fractional nature of its negative argument. Notice also that 
there is a singularity oc (A — l)~(2 + at A = 1. To the left (right) of A = 1 this 



5 In one dimension all bubble diagrams are finite just like the integral over thermal phase space 
in Eq. ([6]) is. The important difference is that the coefficients of the expansion in A are alternating 
in sign [TJl [13] while they are strictly positive in Eq. ([6]) . 
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singularity looks like a real (imaginary), nonintegrable pole which follows from the 
limit formula [TJ] 

lim Li s (e") = r(l - ^(-/O- 1 , (a = - (| + A < 1) . (10) 

3.2 Inverse Borel transformation for A < 

The inverse Borel transform P mass (A) of Pp maas (A) is defined as 

P mass (A) = J>£(X) = / dte-'Bp^JXt) , (11) 

where 

/oo 
dte-*Li_(s +l )(A*). (12) 

The following integral representation of Li s (2;) holds for all complex s and z |16j : 

iz f i {—z) u , s 

Li s z = — / dit- \ ; - — -, 13 

K ' 2 J c (1 + u) s sin(Trit) ' v ; 

where the path C is along the imaginary axis from — ioo to +ioo with an indentation 
to the left of the origin. Inserting Eq. (jl~3l into Eq. (JT2l) for A = — |A| < and 
interchanging the order of integration, we have 

Pi{\) = -i f du S + ^ e -^t. e d-H.)iog(-A) r{u + 2) (14) 

Since, by Stirling's formula[§ the gamma function r(w + 2) decays exponentially fast 
for u — > iz oo , the integral over it in Eq. (TBI) exists and defines the real-analytic 
function P t ( A), (A < 0). 

We have performed the inverse Borel transformation in Eq. (j!2p numerically for 
Z = 0, 1 L 2, 3, 4. In Fig.d we have plotted the functions P(A) , (Z = 0, 1, 2, 3, 4) , for 
—15 < A < 0. Modulo powers in it and for \u\ ^> 1 the integrand in Eq. (I14R roughly 



is of the form 

c cos(log(— X)x) e~ ax + s sin(log(— X)x) e~ ax , (c, a real , a > 0; s imaginary) , (15) 



6 r(z) = V2^z z ~i e~ z e H ^ where H(z) = J2„>o ({z + n+ 1/2) log (l + - 1) converg. 

for z 6 C_ and liniM_ ) . 0o H(z) = 0, see [T7] , 

7 This follows from Eq. (Tl2|) and the fact that Im Li_( 3 + ;) ( z ) 



= for z < 0. 
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Figure 2: The functions P,(A) for I = (a), I = 1 (b), I = 2 (c), I = 3 (d), I = 4 (e) 
and — 15 < A < 0. Amusingly, / measures the (increasing) time at which snapshots 
are taken of an 'ocean wave' crashing against the 'beach' at A = where Pi = 0. 
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-0.1 -0.05 



Figure 3: The behavior of the function RePo(A) hi the vicinity of A = 0. 

where x denotes the positive-real integration variable. Thus a fit function $/(A) of 
the following form is suggested^: 



Ef=0 Q 2 r+l(log(-72r+lA) ) 

Efi ^(log(-^A)) 



2r+l 



(17) 



where the j2r+i, <>2s are positive- real, and «2r+i) fos are real. The reasons for intro- 
ducing higher (even and odd) powers of log(— A) into the numerator and denominator 
of Eq. ffTTl) (not setting 5 2s = 72r+i = 1 and Si — R\ + 1 = 1) are due to the con- 
tributions to the integral in Eq. (TT^j) from moderate values of |u| and the presence 
of the factor (1 + u) 3 ^ 2+l . We have checked that the functions displayed in Fig. [2, 
indeed, fall off logarithmically slowly for — A — > oo. Except for the branch cut along 
the positive-real axis is a meromorphic function of A. As such it can be continued 
arbitrarily close to the cut. 



3.3 The case A > 

The functions P(A) exhibit a branch cut along the positive-real axis. This is because 
of their dependence on powers of log (—A) as suggested by Eq. dHJ. At A = there is 
a cusp in Re P,(A) with Re P(A = 0) = 0, see Fig.EJ The slope of Re F,(A) at A = ^0 
is ± oo. Excellent fits reveal the following expressions for $; , (/ = 0, 1, 2, 3, 4, 5): 



s Onc has: 



dxe ax sin(log(-A):c) = - , 

a 2 + (log(-A)) 2 

dx e— cos(log(-A» = - ^~\. 2 , (a > 0) . 

a 2 + (log(-A)) 2 



(16) 
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0.001 
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(e) 

Figure 4: The real (black) and the imaginary (grey) part of the functions Pi(X) for 
I = (a), I = 1 (b), I = 2 (c), I = 3 (d), I = 4 (e) and positive A. 



* (A) = 0.0570 M-0.1B4X) 



*i(A) 



$ 5 (A) 



$ 3 (A) = -0.0544 
$ 4 (A) = 



1 + 0.220(log(-0.494A)) 2 
0.0212 log(-10.2A) + 0.00142(log(-0.109A)) 3 
1 + 0.128(log(-1.09A)) 2 + 0.0544(log(-0.886A)) 4 
log(-0.367A) 



l + 0.641(log(-0.513A)) 2 ' 
-0.0722 log(-1.66A) + 0.00780(log(-1.95A)) ; 



1 + 0.212(log(-1.90A)) 2 + 0.0864(log(-1.06A)) 4 

14.71og(-0.505A) - 1.01(log(-0.917A)) 3 - 0.0172(log(-0.0568A)) 5 
1 + 5.84(log(-0.0156A)) 2 + 0.925(log(-10.3A)) 4 + 1.68(log(-0.641A)) 6 

(18 



Fig. H] shows for A > plots of the Re$;, which are continuous across the cut, and 
Im$2, whose respective signs are chosen such as to minimize the value of A > 
where they first intersect with Re $j (approach to the cut as A = Re A ± i 0). Notice 
that Re and Im vanish precisely at A = 0. Notice also that for A > and in the 
vicinity of A = the modulus of Re grows much more rapidly than the modulus 
of Im $/. 
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3.4 Interpretation of the results of Sec. 13.31 



Due to the meromorphic nature of the functions $/ away from the cut along the 
positive-real axis and the fact that Re is continuous across this cut we have 
no choice but to regard, up to a modification involving finitely many termg^l, the 
quantity 

P mass (/3) = 3, Re J2 a i p i (KP)) (19) 

as a unique prediction for the pressure exerted by the massive modes in the confining 
phase of an SU(2) Yang-Mills theory, compare with Eq. ([6]). We expect this pressure 
to be positive and rapidly growing for A > 0, see j9] where the asymptotic expansion 
of Eq. (jSJ) was used to predict the temperature dependence of the pressure for an 
'electron' gas. For accurate numerical predictions the integer coefficients a\ would 
have to be known. In a A0 4 -theory the algebraic dependence of the number of 
connected vacuum bubble diagrams on n should be extractable from a numerical 
analysis of the ground-state energy of an anharmonic oscillator, see [T2], [13] . 

It is important to realize that at ft = oo (or T = or A = 0) we have P mass (/3 = 
oo) = 0. That is, the contribution to the cosmological constant arising from an SU(2) 
Yang-Mills theory at zero temperature is nil. Thus the according 'tree-level' result 
[H [2] , which relates to the zero of the effective potential for the vortex condensate 
(oc the 't Hooft-loop expectation [H [2]) survives the integration of fluctuations as 
performed in the present work. 

But what about the imaginary part in P mass ? As is readily observed from Fig.H] 
(and is easily proven analytically for all I) the modulus of Re$/(A) grows much 
faster than |Im$/(A)| with increasing A starting at A = 0. Thus for a certain range 
of small valued (3~ l the pressure P mass (/3) is dominated by the real part. 

The presence of a growing imaginary contamination signals an increasing devia- 
tion of the system from a thermodynamical behavioi0. Namely, we interprete the 
occurrence of a sign-indefinite imaginary part0 as an indication for the violation 
of a basic thermodynamical property: Spatial homogeneity. That is, the growth of 
the ratio R = ImP mass /ReP mass with increasing temperature is a measure for the 
increasing importance of turbulence-like phenomena in the plasma. At R = 1 the 
thermodynamical description of the system fails badly; the system then is highly 
'nervous' and close to the Hagedorn transition. 

9 We have used the large- n expressions for the coefficients N n . Precise predictions would have 
to subtract the first few terms from the expression in Eq. (|19|) (evaluated with large-n coefficients) 
and add terms with realistic multiplicities. 

10 For example, A = 0.1; .01; 0.0004 corresponds to 

= £ = - (log (sly)) = 0.173; 0.124; 0.0886, respectively. 

n No well-defined partition function directly generates an imaginary part for the pressure. 
12 Imaginary contributions to the pressure lead to localized exponential built-up and collapse of 
energy density about the equilibrium situation, that is, to turbulences. 
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We believe that the so-predicted occurrence of sizably nonthermal behavior at 
sufficiently large temperature is the reason for the failure of stabilization of the 
magnetically confined plasma in tokamak experiments. This presumes that the 
description of the electron and its neutrino and their quantum mechanics [18J is based 
on a pure SU(2) Yang-Mills theory of scale A = m e = 511 keV, see [H [21 [9]. Notice 
that the electron and its neutrino are the only stable excitations in the confining 
phase of this theory, compare with Fig.[TJ The occurrence of poorly understood 
microturbulences and internal transport barriers in tokamaks was reported for ? ~ 
0.05 ... 0.09, see JET's results and ITER's design report. 

4 Summary and Outlook 

In the present work we have performed an in-principle calculation of the pressure 
El of an SU(2) Yang-Mills theory at low temperature. In particular, we have shown 
that the pressure vanishes precisely at zero temperature. We have also demonstrated 
that nonthermal effects start to become sizable at a certain critical temperature T*. 
Numerical predictions for T* and the dependence of the pressure on temperature need 
input about the precise algebraic dependence of the number of connected vacuum 
bubbles on the number n ^> 1 of vertices in a A(/> 4 -theory. Also, one would need 
to know these number^* 4 ! for the first few n. We hope that this information will be 
available soon. 

The tool employed to derive the above results is to show Borel summability of 
an asymptotic-series representation of the pressure at negative values of a coupling 
constant A suitably defined for T > 0. The so-obtained functional dependence 
is meromorphic in the entire complex A-plane except for a branch cut along the 
positive-real axis. The crucial observation is that the real part of this function 
is continuous across the cut and that its modulus grows faster than that of the 
imaginary part for a certain range of increasing A > 0. The presence of the imagi- 
nary contamination signals a deviation from thermal equilibrium by local violations 
of spatial homogeneity (plasma turbulences). We believe that our results are of 
relevance in addressing the observed but poorly understood microturbulences and 
internal transport barriers ocurring in terrestial fusion experiments with magnetic 
plasma confinement. 

The generalization of our SU(2) results to SU(3) is trivial since for SU(3) the 
number of vortex-loop species at a given time is just twice that of the SU(2) case. 

13 Other thermodynamical quantities are derivable from the real part of the pressure by Legendre 
transformations. 

14 They are known for n < 6 [14] . 
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